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INTRODUCTION 

The Lageos satellite is in a high-altitude (5900 km), almost circular orbit about the earth. 
The orbit is retrograde: the orbital plane is tipped by about 1 10 degrees to the earth’s equatorial 
plane. The satellite itself consists of two aluminum hemispheres bolted to a cylindrical beryllium 
copper core. Its outer surface is studded with laser retroreflectors. For more information about 
Lageos and its orbit see Smith and Dunn (1980), Johnson et al. (1976), and the Lageos special 
issue (Journal of Geophysical Research, 90, Bll, September 30, 1985). For a photograph see 
Rubincam and Weiss (1986) and a structural drawing see Cohen and Smith (1985). Note that the 
core is beryllium copper (Johnson et al., 1976), and not brass as stated by Cohen and Smith 
(1985) and Rubincam (1982). See Table 1 of this paper for other parameters relevant to Lageos 
and the study presented here. 

Lageos’ orbit is the most accurately modeled of any satellite (Cohen and Smith, 1985). 
However, after subtracting out most of the known forces acting on the satellite, such as the gravi- 
tational attraction of the sun and moon, direct solar radiation pressure, etc., there is still a re- 
sidual along-track acceleration which remains to be explained. That acceleration is the subject of 
this paper. 

The residual is given as a function of time in Figure 1 . The figure shows the average 
monthly values of S, the unexplained along-track acceleration. It clearly acts like drag and has 
a mean value of -3.33 x 1 O' 12 m s' 2 . It brings Lageos closer to earth by 1 .2 mm per day. 
Moreover, there are fluctuations in the acceleration which can be as large as the mean value. At 
times S drops almost to zero, as in March of 1983. Most of the fluctuations are obviously corre- 
lated with the sun-orbit geometry: the largest ones occur when Lageos spends time in the earth’s 
shadow. 

The problem is to account for the curve shown in Figure 1 . Afonso et al. (1980), Mignard 
(1981), Rubincam (1982), Afonso et al. (1985), and Earlier et al. (1986) all invoke charged par- 
ticle drag to explain the average decay. Neutral atmospheric drag appears to explain only about 
10 percent (Rubincam, 1982; Afonso et al., 1985). Anselmo et al. (1983) and Barber et al., 
(1986) invoke earth-reflected sunlight to account for the fluctuations. The reflecting surface can- 
not be Lambertian (Barlier et al., 1986; Rubincam and Weiss, 1986). Morgan (1984) invokes 
earth-reflected sunlight to explain both the average drag and the fluctuations. 
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The purpose of this paper is to evaluate a new mechanism to explain the observed average 
drag: thermal drag. It is probably best thought of as a variant of the Yarkovsky effect (Opik, 
1951; Bums et al., 1979; and Rubincam, 1982). 

The basic idea is this: the earth emits infrared radiation which is ultimately due to solar heat- 
ing. The exitance is a uniform 232 W m' 2 over the earth’s surface, to a first approximation 
(Stephens et al., 1981, Table 4a; Sehnal, 1981, p. 169). Naturally, Lageos will intercept some of 
this radiation. 

Assume Lageos’ spin axis is in the plane of the orbit, as shown in Figure 2; the orbit is taken 
to be circular. Because Lageos’ spin axis is fixed in inertial space, first the retroreflectors on one 
hemisphere and then the other will be alternately heated by the earth’s infrared radiation as the 
satellite circles the earth. (So will the aluminum surface, but this will be ignored here.) The alter- 
nate heating of hemispheres causes a temperature asymmetry between them. The hotter hemis- 
phere radiates away more energy and thus more momentum than the cooler. The satellite there- 
fore feels a force along the spin axis directed away from the hotter hemisphere, due to momem- 
tum conservation. 

If the retroreflectors had no thermal inertia, then a hemisphere would be hottest when the 
earth is directly over its pole. In this case the net along-track component of the force is zero be- 
cause of cancellations when averaging over one revolution. However, since the retroreflectors do 
have thermal inertia, a hemisphere becomes hottest after the earth has passed over its pole; there 
is a delay. It is then easy to show that the force does have a net along-track component when av- 
eraging over one revolution (and is best envisioned when the lag is 90 degrees; see Figure 2). 
Moreover, the average along-track acceleration is always opposite to the direction of motion; in 
other words, it acts like drag and causes the orbit to decay. 

This paper estimates the magnitude of the along-track acceleration due to the thermal drag 
by modeling the thermal behavior of the retroreflectors. The result is that the thermal drag 
accounts for about 47 per cent of the observed average drag. A rough estimate of the aluminum 
surface’s contribution is put at less than 5 per cent. In view of these results, thermal drag may be 
the dominant drag mechanism operating on Lageos; the other proposed mechanisms mentioned 
above may have to be reassessed. However, thermal drag apparently cannot explain the large 
fluctuations in the observed along-track acceleration. Some other mechanism must be invoked to 
account for them, such as earth-reflected sunlight or perhaps some unthought-of phenomemon. 
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entry and exit: “S-N”, for instance, means Lageos entered the shadow from the south and exited from the north. The 
dashed line gives the average acceleration: -3.33 x 10' 12 m s~ 2 . 
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Figure 2. Schematic diagram of the infrared effect when the spin axis is in the plane of the orbit. 
The earth’s infrared radiation (wavy arrows) heat up the retroreflectors (not shown) on Lageos’ 
surface, causing a net force (thick arrows) along the direction of the spin axis. The thermal lag 
angle is taken to be 90 degrees in this diagram merely to illustrate that there is an along-track 
deceleration. The actual angle is closer to 37 degrees. 


THE RETROREFLECTORS 


There are 426 comer-cube laser retroreflectors on Lageos. They cover 42 per cent of Lageos’ surface 
(Johnson et al. , 1976). Of these 422 are made of fused silica. The remaining four are made of germanium 
(Johnson et al., 1976). Only the silica retroreflectors are discussed here. The few made of germanium are not 
considered at all in this paper. 

The relevant properties of a silica retroreflector are given in Table 1 . The values for the density p, specific 
heat C p , thermal conductivity k, infrared emissivity e, and average temperature T„ come from Bendix (1974, ap- 
pendices A and K). The diameter of the outer face d comes from Johnson et al. (1976), while the height h from 
the face to the comer is estimated from their Figure 5, which shows a drawing of a retroreflector. The effective 
radius R r of a sphere with about the same volume as as retroreflector comes from approximating the shape as a 
cone with diameter d and height h: ir d 2 h/12 = 4ir R r 3 /3. The fraction f c of the face area to the total area comes 
from measurements made on a cardboard model. The values for R r and f c are thus somewhat crude, but should 
be good enough for the purposes considered here. 


RETROREFLECTOR MODEL 

The first step in finding the along-track acceleration S from the thermal drag is to thermally model a retrore- 
flector. This is done as follows: the retroreflector is assumed to be a sphere with radius R r and illuminated by 
three sources: infrared radiation from the earth; sunlight; and infrared radiation from the aluminum cavity in 
which the retroreflector sits. The retroreflector is taken to be conductively isolated from the aluminum (Johnson 


Table 1 

Constants relating to Lageos, the earth, and the universe. 


Symbol 

Quantity 

Numerical value 

a 

orbital semimajor axis 

1.227 x 10 7 m 

Cp 

retroreflector specific heat 

712 Jkg 1 K 1 

fo 

ratio of face to total area of retroreflector 

0.33 

I 

orbital inclination to earth’s equator 

109.9 degrees 

M L 

mass of Lageos 

407 kg 

n L 

Lageos mean motion 

4.65 x 10‘ 4 s’ 1 

Rl 

radius of Lageos 

0.3 m 

Rr 

retroreflector radius 

0.01331 m 

T 0 

retroreflector temperature 

276.7 K 

€ 

retroreflector emissivity 

0.9 

K 

retroreflector thermal conductivity 

1.34 WK' 1 m 1 

P 

retroreflector density 

2200 kg m" 3 

af ir ° 

sub-earth irradiance at Lageos altitude 

62.55 Wm' 2 

Re 

radius of earth 

6.371 x 10 6 m 

c 

speed of light 

2.9979 x 10 8 m s’ 2 

a 

Stefan-Boltzmann constant 

5.67 x 10' 8 W m’ 2 K" 4 
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in the above equation gives 


AT(R r ) 


AF IR f 0 

4oT 0 3 + 1 v P C p R R 
3e 


or 


At (R r ) = 


AF ir f 0 e' iS 
4aT 0 3 (1 + Z 2 )' /2 


( 10 ) 


where 


8 = Arc tan £ 


and 


£ _ V P Cp^R 
12 ecxT 0 3 


( 11 ) 


( 12 ) 


The approximations for sin w and cos w will be justified later. Interestingly, k drops out of (10). Also, note that 
8 is positive or negative depending on whether v is positive or negative. 


TEMPERATURE DIFFERENCES BETWEEN RETROREFLECTORS 

Equation (10) gives the expression for the surface temperature for a given retroreflector. The next step is to 
find how AT(R r ) varies from retroreflector to retroreflector due to differences in their illumination by the earth. 

Assume for simplicity that the earth is infinitely far away from Lageos; then the impinging rays come in 
parallel lines and 


AF ir — AF ir ° D (0, X, 0 E , X E ) 
where AF IR ° is the irradiance at the sub-earth point and 

oo + L , 

D(0,X,0 E ,X E ) = 4tt 2 2 ( ~ W* OeAe) Y l j (0,X) 

(Rubincam and Weiss, 1986). Here the Y L J (0,X) are the spherical harmonics used in quantum mechanics; the 
asterisk (*) means complex conjugate. The coordinates (0,X) are the colatitude and longitude, respectively, of a 
retroreflector on the satellite (see Figure 3). The angle 0 is measured from the spin axis, which is the satellite z L 
axis. The coordinates of the earth are (0 E ,X E ) in this system. The d L are coefficients given in Table IV of 
Rubincam and Weiss (1986). AF IR ° is easy to find; due to inverse square attenuation AF IR ° = 232 (R E /a) 2 = 
62.55 W m" 2 , where R E is the earth’s radius and a is the orbital semimajor axis. 

Consider only the L - 1 , J =0 term in D(0,X,0 E ,X E ); then d, = Vi and D (0,X,0 E , X E ) = cos0 cos0 E /2. 
The change in temperature thus depends only on colatitude 0 on the satellite. The dependence on longitude X is 
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Figure 3. Coordinate system (x L , y L , z L ) rigidly attached to Lageos. The z L axis points along the spin axis of the 
satellite. The retroreflector is located at (9, A.), whre 9 is colatitude and A is longitude. The earth is at (0 E , A E ). 


ignored and assumed to be small due to the rapid rotation of Lageos (a spin period of about 7 s at the end of 1983; 
E. M. Gaposchkin, private communication, 1985). The temperature also depends on time through 0 E . In fact, the 
time dependence can be made explicit by writing 


cos 0 E = 


e in L t + e -in L t 


2 


0E — n L* 


without loss of generality. Here n L is the mean motion of the satellite and its value is given in Table 1. 
Using these data in (10) - (12) give 

AT = B cos 0 cos (n L t - 8) 


(13) 


1 
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where n L takes the place of v, and where 


AF [R °f 0 

B = = 1.72 K (14) 

8aT 0 3 (l + £ 2 )' /2 

and 

5 = Arc tan (0.7465) = 36.7 degrees (15) 

The numerical values come from using Table 1 . The cos0 factor in (13) shows that there is an axial asymmetry in 
the retroreflector temperatures across Lageos; one hemisphere is hotter than the other. The cos (n L t - 8) in the 
same equation shows that the hemispheres alternate in hotness. Equations (13) - (15) thus give the quantitative 
formulation of the qualitative statements made in the Introduction and are fundamental to this paper. Note that by 
(8) (kR R ) 2 = -0.096i, justifying the approximation for sin w and cos w made earlier. Note also that the lag angle 
given by (15) corresponds to about 23 minutes of time. This is long compared to the spin period and justifies ig- 
noring the temperature dependence on It is short compared to the orbital period of 3.76 hours. 

MAXIMUM DRAG 

The last step is to find the along-track acceleration S due to the temperature asymmetry. Consider the flat 
outer surface of an actual retroreflector, so that the assumption of sphericity is dropped; equations (13) - (15) 
will, however, still be assumed to be valid. It is easy to show that if the face has area dA and temperature T and 
emits radiation according to Lambert’s law (e.g., Brown, 1965, p. 225), then the photons leaving it carry away 
momentum at the rate of 


^ = 2eoT 4 dA/3c 

each second normal to the surface; c is the speed of light. Now because B << T 0 , by (3) and (4) the component 
of this force in the satellite z L direction (along the spin axis) is approximately 

(T c 4 + 4T 0 3 AT) cos 6 dA 

dt 3c 

The T 0 4 term drops out when integrating over Lageos’ surface, leaving 


... 4ireR L 2 AF IR °f 0 

z L = - (0.41) cos (n L t-8) 


9M l c (1 + £ ) 


2x1/2 


(16) 


as the total acceleration due to all of the retroreflectors. The factor 0.41 comes from the fact that the silica 
retroreflectors cover only 41 percent of Lageos’ surface. 

The along-track acceleration S is z L sin 0 E . Averaging S over one revolution gives 


<S> 


max 


(0.41) 


2tt 6R l 2 AF 1R °f 0 
9 M l c ( 1 + £ 2 ) 1/2 


sin 8 


(17) 
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for the case of the spin axis lying in the orbital plane. Substituting numerical values in this equation give 

<S>max = - 1.877 X 10 12 m s' 2 ( 18 ) 

This is the extreme value <S> achieves for the thermal model considered here. Any other spin-orbit geometry 
will give a smaller drag. 

The value given by (18) is about 56 per cent the observed average drag. Hence the thermal drag becomes a 
major contender for explaining the average acceleration as shown in Figure 1 provided the spin-orbit geometry is 
favorable. The subject of geometry is taken up next. 


LAGEOS ALONG-TRACK ACCELERATION 

What is desired now is a general expression for <S> when the spin axis is not necessarily in the plane of the 
orbit. The purpose is to see how <S> varies over long periods of time as the spin axis shifts and the orbit 
precesses. 

Consider the coordinate system shown in Figure 4. The origin is at the center of the earth and the z-axis 
pierces the North Pole. The x-axis points to the vernal equinox. Assume once again that Lageos’ orbit is circular; 
then Lageos will have a unit position vector (Goldstein, 1950), p. 109.): 

r = [cosftcos(to + f) - cos I sinftsin(to + f)] x 

+ [sinftcos(to + f) -I- cos I cosftsin(to + f)] y (19) 

+ [sin I sin (to + f)] z 

in this system. Here I is the inclination of the orbit to the earth’s equator (a constant 109.9 degrees), ft the nodal 
position in the equatorial plane, to the argument of perigee, and f the true anomaly. The unit vectors x, y, z, lie 
along their respective axes. The unit along-track vector t is 

t = [-cosftsin(w -I- f) - cos I sinftcos(o) + f>] x 

+ [-sinftsin(<o + f) + cos I cosft cos(o> + f)] y 

+ [sin I cos(o) + f)] z 

The unit vector in the spin direction is 


s = s x x + s y y + s z z 


where 


S x 2 + S y 2 + s z 2 = 1 

The acceleration will lie along the s-axis and be proportional to cosi|j, i.e. r « cos i|i s , where cos 4* = r d • s. 
Here r d is the delayed unit position vector in which to + f - 8 replaces to + f in (19) to allow for the thermal lag. 


11 




The along-track acceleration S is then obtained from r via S = 7* t. Using (19) and the succeeding equations and 
averaging co + f from 0 to 2 tt give 


<S> = <S> max [l-s z 2 + Vi (3s z 2 -1) sin 2 1 

+ s z sin21 (s y cosft - s x sinft) 

( 20 ) 

+ Vz (s x 2 - s y 2 ) sin 2 I cos 2ft 
+ s x s y sin 2 I sin 2ft] 


The time dependence of (20) comes in implicitly through s x s y ,s z , and ft. These quantities change only 
slowly compared to 10 + f and were held fixed in the averaging process. The explicit dependence of ft on t can 
be written 

ft = ft (t-t 0 ) + ft 0 


where 


ft = 0.3425 degrees/day 

ft 0 = 28.5596 degrees 

t 0 = 42902.5 (Modified Julian Date) 

so that the nodal position moves at a constant rate along the equator. The dependence of s x , s y , and s z on t is a 
more complicated story. One must find their initial values at launch (t = t D ) and then find how the torques on the 
satellite move the spin axis. 

The initial values of s x , s y , and s z can be found from the flight data (J. D. Kraft, oral communication, 1978). 
Lageos was launched on 4 May 1976 at 4 AM Eastern Daylight Time. Apogee kick motor separation occurred 
5457 seconds later at +4.5 degrees geocentric latitude and +20.8 degrees geocentric longitude. The elevation of 
the spin axis was -1 1 .0 degrees and the azimuth was + 158.7 degrees in the topocentric frame. These data yield 

s x = ■+■ 0.25433 

s y = ± 0.27650 on 4 May 1976 
s z = + 0.92675 

for the initial values in the frame of Figure 4. The details of the calculations, which involve the multiplication of 
rotation matrices, are omitted here. The ± signs appear because of uncertainty in the sense of rotation of the 
spin. They create no problems in evaluating (20), since that expression is invariant under change of sign of s. 

This merely reaffirms that it is the direction of the spin axis which is important, not the spin vector. 
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Substituting the above values in (20) gives 


<S> = <S> max . [ + 0.8381 

- 0.2229 cos (ft - 222.6 deg) (21) 

- 0.0624 cos (2ft - 85.2 deg)] 

This expression shows that the value for <S> given by (18) must be reduced by a factor of 0.8381: 

<S>sec = -1 -57 x 10' 12 m s' 2 (22) 

to obtain the secular part of S immediately after launch. This is 47 per cent of the observed average of -3.33 x 10" 12 
m s' 2 . Also, <S> varies sinusoidally with frequencies ft and 2ft, with the cosft term being the principal one. 

Since the factor in brackets in (21) becomes as small as 0.8381-0.2229-0.0624 = 0.5528, <S> varies from (22) 
by 100 (0. 8381-0. 5528)/0. 8381 = 34 per cent. 

The spin axis will change its position in space as Lageos is despun by the earth’s magnetic field. The gravi- 
tational torque on Lageos is negligible (Barlier et al., 1986). The Lageos spin axis should eventually line up with 
the earth’s spin axis, since the average magnetic field points north-south. In this case s x = s y = 0 and s z = 1 in 
(20), giving 


<S> = -1.66 x 10' 12 m s' 2 

with no periodic terms. This is not much different from (22); hence the average drag will probably not change 
much as the Lageos spin axis aligns itself with the earth’s spin axis. However, no detailed treatment of this 
problem is given here; both the precession of Lageos’ spin axis and the time scale associated with despin have 
been ignored, as well as the evolution of the orbit when Lageos spins very slowly. These will certainly be of fu- 
ture interest. 


DISCUSSION 

The thermal model presented here explains a little less than half the observed average drag. The question 
naturally arises as to whether the thermal drag explains nearly all of it, or whether a better model will give an 
even smaller contribution to the drag than that computed here. This obviously cannot be answered without study- 
ing a more detailed model; so it is worthwhile to discuss in what areas further progress might be made towards a 
new model . 

One of the simplifying assumptions made here was that all of the retroreflectors have the same temperature 
T 0 in the absence of heating by the earth. But a worst case analysis where the Lageos spin axis points towards the 
sun gives a temperature difference of 24 K between the antipodal retroreflectors (Bendix 1974, Appendix K). 
Should this be of concern in a new model? 

The answer is probably no, for two reasons. One is that for the initial spin axis orientation the sun never de- 
parts from Lageos’ equatorial plane by more than 15 degrees (details of the computation are omitted); and never 
more than 23.5 degrees when the spin axis aligns itself with the earth’s axis. Hence both hemispheres are about 
equally illuminated — far from the worst case where one is in shadow and the other in sunlight. So the differ- 
ences between retroreflector temperatures are probably only a few degrees. The other reason is that <S> as 
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given by (17) is relatively insensitive to changes in T 0 when £ is near 1 . Thus it appears that there is not much to 
be gained by looking at a temporally- and spatially-varying T 0 . 

In assuming the retroreflectors were spherical only the / = 0 term in (7) was considered and the higher har- 
monics were ignored. Carrying through the analysis for / = 1 shows a negligible effect on <S>; so the neglect 
of the l > 0 terms does not appear to be serious. 

Also ignored here is sunlight reflected from the earth; this too will heat up the retroreflectors. But the effect 
is tiny. Going from no sunlight to full sunlight changes a retroreflector’s temperature by only ~ 3 K for a cavity 
temperature of 303 K (NASA, 1975, Table 3-1). (Lageos’ actual cavity temperature is about 328 K; see Bendix, 
1974, Appendix K.) The irradiance due to reflected sunlight is about (0.3) (1376)/(4*4) ~ 26 W m 2 at Lageos’ 
altitude. Here 0.3 is the earth’s albedo (Stephens et al., 1981) and 1376 W m‘ 2 is the irradiance from the sun 
(Hickey et al., 1980). One factor of 4 comes from spreading the incident sunlight over the earth’s surface, while 
the other comes from the inverse- square law. Therefore the temperature change due to reflected sunlight must be 
on the order of (26) (3)/l 376 = 0.06 K, which is small. 

The earth was assumed here to be a point source illuminating Lageos. Actually, the earth subtends a diame- 
ter of about 60 degrees as viewed from Lageos, so that more than half its surface is heated by the earth. The extra 
illumination is offset somewhat by specular reflection (ignored here) off the retroreflectors with large zenith 
angles, but the matter may be worth further consideration. 

Probably the most progress will be made in giving a retroreflector a more realistic shape than a sphere and 
modeling its interaction (via radiation) with the cavity. NASA (1975) describes a numerical model consisting of 
66 nodes where the retroreflector is shaped like a cone and sits in a cylindrical cavity. This perhaps should be the 
point of departure for a future model, and it will be interesting to see whether the front face temperature changes 
go up or down; up means higher drag. 

There is some evidence that they will go up. Experimental tests show that infrared radiation increases the 
steady-state front face temperature for the sub-earth retroreflector by about 10 K, for a cavity temperature of 303 
K (NASA, 1975, Table 3-1). The time-variable case will divide this amplitude by 2(1 + £ 2 ) 1/2 , the factor of 2 
coming from D(6, \, 0 E , X. E ), and the (1 + £ 2 ) l/2 coming from (10). This would give a value of B in (14) of 
~4K, assuming the same lag angle, which would about double the drag and explain all of the observed secular 
decrease in the semimajor axis of the orbit. This clearly warrants further investigation. 

Nothing has been said so far about the contribution by the aluminum shell of Lageos to the thermal drag. 

Its contribution can be shown to be small, as evidenced by the following argument. The worse case analysis by 
Bendix (1974) gives a maximum temperature difference of 8 K between the two aluminum hemispheres for solar 
heating. This would make the amplitude 4 K instead of 1 .72 K in (13) for sunlight. Since the earth’s irradiance 
(62.55 W m' 2 ) is 22 times smaller than the solar irradiance (1376 W m‘ 2 ; Hickey et al. , 1980), the amplitude for 
the earth would be approximately 4/22 = 0.18 K, which is about 10 per cent that of the retroreflectors, assuming 
the same lag angle. The actual case must be even smaller than the worst case. Hence Lageos’ aluminum surface 
explains only about 5 per cent or less of the observed average drag and plays only a small role in the thermal 
drag. 

The model described here cannot explain the observed fluctuations shown in Figure 1 . The model gives 
sinusoidal fluctuations in drag with a period of 360/0 = 1051 days (see equation 21). Its signature is entirely dif- 
ferent from the observed one: the biggest observed fluctuations clearly occur when the orbit intersects the earth’s 
shadow. Could diurnal variations in the earth’s exitance somehow account for them? 

The answer to this question appears to be no. Clouds and ocean maintain a fairly constant exitance through- 
out the day, according to METEOSAT data (Saunders and Hunt, 1980). The Sahara seems to give the most 


15 


extreme variation in exitance; but this still only amounts to a change of about 50 Wm" 2 from day to night (Figure 
2 of Saunders and Hunt, 1980). This will be reduced by at least a factor of 4 due to inverse- square or higher at- 
tenuation, making it small compared to the uniform 62.55 W m" 2 . Averaging over many rotations of the earth 
will make it even smaller, since Lageos is not in an orbit resonant with the earth’s rotation. Explanations for the 
observed fluctuations must be looked for elsewhere. Anselmo et al. (1983) and Barlier et al. (1986) propose radi- 
ation pressure from sunlight reflected off of the earth. 

The thermal drag depends crucially on the spin-orbit geometry. If the spin axis is normal to the orbital 
plane, then there is no drag. If the spin axis is in the plane of the orbit as in Figure 2, then there is maximum 
drag. These facts allow predictions to be made. If the thermal drag is the main cause of the secular decay of 
Lageos’ orbit, then Doppler radar observations should confirm that Lageos’ spin axis is never perpendicular to 
the orbital plane when sizable drag occurs. Also, the drag on Lageos II should be calculable if the spin axis orie- 
nation is known. (Lageos II should differ from Lageos only in having an orbital inclination of 59 degrees instead 
of 1 10 degreees, assuming a normal launch. Its launch date is now uncertain, due to the Challenger disaster.) 

If the thermal drag is the dominant drag mechanism operating on Lageos, then charged particle drag must 
be smaller than previously believed (Afonso et al., 1980; Mignard, 1980; Rubincam, 1982; Afonso et al., 1985; 
and Barlier et al., 1986). In any case, the heating of Lageos’ retroreflectors by the earth’s infrared radiation ap- 
pears to play a significant role in the secular decrease in the semimajor axis of Lageos’ orbit. 
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